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Abstract: We derive sum rules which constrain the spectral density correspond- 
ing to the retarded propagator of the T^y component of the stress tensor for three 
gravitational duals. The shear sum rule is obtained for the gravitational dual of the 
jV = 4 Yang-Mills, theory of the M2-branes and M5-branes all at finite chemical 
potential. We show that at finite chemical potential there are additional terms in 
the sum rule which involve the chemical potential. These modifications are shown 
to be due to the presence of scalars in the operator product expansion of the stress 
tensor which have non-trivial vacuum expectation values at finite chemical potential. 
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1. Introduction 

Sum rules play an important role in constraining spectral densities of strongly coupled 
fluids. There are important sum rules which constrain spectral densities of stress 
tensor in QCD |l|, ^, § and the Ferrell-Glover-Tinkham sum rule satisfled 
by current-current correlation function in the BCS superconductor P, 0]. The gauge 
gravity duality provides a framework for evaluating the two point functions of various 
conserved currents of strongly coupled theories which admit a gravity dual. Thus one 
can use this framework to obtain sum rules. Obtaining sum rules within this frame 
work provides information of the analytic structure of the Green's functions in the 
strongly coupled limit. Romatschke and Son derived two sum rules for the stress 
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tensor two point function for strongly coupled A/" = 4 Yang-Mills theory. This was 
then generalized for Chamblin-Reall backgrounds which are dual to non-conformal 
theories in |T0|. The sum rule for the R-charge correlator in strongly coupled 



TV = 4 Yang-Mills was obtained in |]T|. 

As emphasized in sum rules are the consequences of analyticity of the 



Green's function in the complex frequency plane. From the field theory point of 
view this results from the unitarity and causality of the field theory. Deriving sum 
rules from gravity provides insight into how the properties of unitarity and causality 
of the boundary field theory are encoded in the gravitational theory. In the 
differential equations which determine the retarded Green's function of interest from 
gravity were studied. The properties of these differential equations were used to 
obtain proof and provide a unified framework to obtain various sum rules in gravity. 

In this paper we obtain the sum rules from gravity for the spectral density 
corresponding to the retarded propagator of the Txy component of the stress tensor. 
This is done for the case of A/" = 4 Yang-Mills, the M2-brane and the M5-brane 
theory in the presence of chemical potential. One of our motivations to examine 
these situations is to determine how the sum rules are modified when the dual theory 
cannot be truncated to pure gravity. The presence of chemical potentials in these 
systems results in additional scalars in the gravitational theory. We will now briefly 
state the result for the A/" = 4 Yang-Mills case. The shear sum rule for Af = 4 
Yang- Mills derived in H in the absence of chemical potential is given by 



r 



-I ^{p{u;)-pT=o{u;)), (1.1) 



where 

p = ImGRiu), (1.2) 

and Gr is the retarded propagator of the T^y component of the stress tensor. pt=o{uj) 
is the spectral density at zero temperature and e is the finite temperature energy 
density. In this paper we examine this sum rule in the presence of chemical potential. 
We find that the sum rule is modified to 



2 N^-K^T^ 
5^ 120 



1 fOO 1 

{{h - k,f + {h - + {h - hf] = - — {p{uj) - PT=,{io)) . 

^ J-oo ^ 

(1.3) 

Where k[s are functions of the three chemical potentials in A^ = 4 Yang-Mills given 
in (|4.6| ). The relation between Tq and the temperature of the Yang-Mills is given 
by ([4.6|) . Note that for the situation when all the charges or chemical potentials are 
equal then the correction vanishes. For this case the the gravity background has no 
additional scalars present. We then show that these additional terms in the sum 
rule are due to the fact that, the operator product expansion of two stress tensor 
involves operators of dimension 4 in addition to the stress tensor. We show that 
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expectation values of these operators precisely have the same dependence in terms 
of the chemical potential to account for the additional terms. 

For the M2-brane theory, the shear sum rule at finite chemical potential is given 

by 

3 v^TT^ N^/'^T^ 

-e H — -{ki + k2-h - k4){ki -k2 + kz - k4){ki - k2 - ks + k^j 

o ziu 

= -/ ^(p(a;)-p^=oM). (1.4) 

Note that for the case of the M2-brane theory, there are 4-chemical potentials which 
are related to the fcj's by ( |5.9|) and Tq is related to the temperature by ( ^.8|) . p is the 
imaginary part of the retarded Green's function shifted by a term which is linear in 
frequency defined in ( |5.24 ). The term proportional to the energy density arises from 



the expectation value of the stress tensor in the OPE, while the term proportional to 
the charge density arises from expectation value of scalars. Finally for the M5-brane 
theory, the shear sum rule at chemical potential remains unchanged in the presence 
of chemical potential. It is given by 

^ (p(^) _ py^p(^)) . (1.5) 

In this case we show the scalars which are turned on due to the chemical potential 
do not have the appropriate conformal dimensions to occur in the OPE of the stress 
tensor and modify the sum rule. 

The organization of the paper is as follows. In the next section we discuss the 
general considerations which go into deriving sum rules in quantum field theories. In 
section 3. we review the proof of the sum rule given in [|lT| for the case of A/" = 4 



Yang-Mills in the absence of chemical potentials. Here we develop a method to 
obtain the behaviour of the retarded Green's function at large frequencies using the 
Fefferman-Graham coordinates. This method also allows us to extract the LHS of the 
sum rule easily. In section 4. we derive the sum rule for D3-branes at finite chemical 
potential and explain the occurrence of additional terms in the sum rule ( p..3|) due 
to the presence of expectation values of scalars. In section 5. we derive the sum rule 
for the M2 and M5-brane theories. Section 6. contains our conclusions. Appendix A 
contains the details of the Fefferman-Graham coordinates which are used to derive 
the sum rule. 

2. Sum rule generalities 

Consider the retarded Greens function corresponding to an operator C of a quantum 
field theory in 4 dimensions given by 

Gnit,x)^temO{t,x),O{0,0)]), (2.1) 
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we then take its Fourier transform given by 

GR{oo,k)= [ d''xe'^'''-''^^GR{t,x). (2.2) 



The spectral density corresponding to this retarded correlator is defined as 

p{oo, k) = j^ {Gr{uj, k) - Gr{uj, k)*} . (2.3) 

It is easy to see that from ( |2.lD and from ( ^.21 ) we see that for bosonic Hermitian 
operators we get 

Gniuj.ky = GR{-uj,-k). (2.4) 
In this paper we will restrict out attention to the case = and we define Gr{uj) = 



Gr{u, 0). From the reality property (p.4|) , we see that in the Taylor series expansion 
of Gr{u) even powers of u have real coefficients and odd powers of u have purely 
imaginary coefficients. Sum rules for the spectral density are obtained by assuming 
certain analytic properties of the Green's function in the complex w-plane. 
The properties are the following 

1. Gr{(jj) is holomorphic in the upper half plane, including the real axis. 

2. lim|i^|_^oo Gr{uj) = if Imw > 0. 

These properties will be referred to as property 1 and property 2 in the rest of the 
paper. From the field theory point of view causality ensures property 1 is satisfied. 
Using these properties, we now indicate the arguments which go into deriving the 
sum rule. By the first property and by Cauchy's theorem we have 

^ / • \ 1 / GR(z)dz . ^, 

Gr{uj + ie) = — i V , (2.5) 

zTii J (J z — oj — le 

^ _ 1 / GR{z)dz 



27Ti jQZ — uj + ie 

for u,e & R and e > 0. The contour G is chosen such that it runs from — r to r for 
some large r and then along the semi-circle in the upper half plane and back to — r 
with r > 0. Because of the second property the integrals along the contour in (|2.5| ) 
reduce to 

Gr{u + ^e) = ^lim f (2.6) 



= — : lim 



GR{z)dz 



27Ti r~^oo z — uj + ie 
Then taking the difference of the two equations in (P^ we obtain 



Gh(.) = lim r !?£^^. (2.7) 

e-s>0+ _^ TT Z — U — le 
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where we have used the definition of the spectral density. Note that the integral runs 
over the real line. The usual sum rule is obtained by evaluating the left hand side of 
( |2.7| ) at w = 0. Thus we have 

GM = hm r (2.8) 

^^o+J_^ n uj-ie 

In many of the situations the Green's functions may not satisfy property 2, that is 
\im\^\^oo G ji{uj) = if Imw > 0. In such situations we subtract divergences and 
obtain a regulated Green's function so that we can still ensure that property 2 will 
be satisfied. The two point function of the stress tensor which we will study in this 
paper falls into this case. 



3. The shear sum rule for uncharged D3-branes 

In this section we review the derivation of the shear sum rule from gravity for = 
4 Yang-Mills at zero chemical potential. This will help us set up notations and 
conventions. To obtain the retarded Green's function of the T^y component of the 
stress tensor one needs to solve the equation of the minimally coupled scalar field 



in the background of the D3-branes [O, 13 . Examining this differential equation 



it is possible to argue that property 1 is satisfied. To ensure property 2, we will 
develop a method to obtain the behaviour of the Green's function in gravity at large 
frequencies. This will enable us to regulate the Green's function and ensure property 
2, which will lead us to the derivation of the sum rule. 

3.1 Green's function from gravity 

As we have mentioned earlier the retarded Green's function of interest in this paper 
is given by 

Gnit^x) = te{t){[T^yit,x),T,yiO,0)]), (3.1) 

where T^y is the xy component of the stress tensor. The dual geometry corresponding 
to A/" = 4 Yang-Mills at finite temperature is the non-extremal D3-brane background. 
The metric of this background is given by 

2 7" 2 

ds' = ^{-fdf + dx" + dy^ + dz^) + —dr\ (3.2) 
The above metric is a solution of the action 
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It is useful to recall the thermodynamic variables of this system. The temperature, 
energy density, pressure and entropy density of this system are given by 



2 



To evaluate the Green's function given in from the gravity dual we first consider 
the metric fluctuation which is dual to the stress tensor T^y. This is given by 

2 

59.y = 0(r)e-^-*+^^^^. (3.5) 

This perturbation obeys the equation of motion of a minimally coupled massless 
scalar in the background (|3.2|) . The A; = mode satisfies the following equation 

+(^ + -] dr<P +^0 = 0, (3.6) 



F2 



where 



(3-7) 



The procedure to obtain the Green's function |]T4| is to first impose in going boundary 
conditions at the horizon r+ and obtain at the boundary r — t- oo. Once this is 
done, the Green's function at temperature T is given by 

contact 

{T) + G 

counter 

{u},T), (3.8) 

Gr{uj) = lim 



87r2L6 



r—^oo 



where Gcontact(7') is the contribution from the contact terms obtained from the on 
shell effective action. The on shell effective action is given by |jTl . 



5„s = P/rf^x(l 



(3.9) 



2' 

where P is the pressure of the solution. Evaluating Gcontact we obtain 

^contact = —P- (3.10) 

G counter {^,T) is the Contribution from the counter terms required to cancel the 
and log(r) divergences in Gii{u). The details of these terms are given in but as 
we will see later, we will not need the detailed structure of these terms. An important 
point to note is that the on-shell action is independent of frequency which ensures 
that the contact term Gcontact is also independent of frequency. 
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Thus the important properties of the Green's function is essentially contained in 
the function 

gR{u) = lim — — . (3.11) 

r— >oo (p 

Therefore to study the behaviour of the retarded Green's function, it is sufficient 
to study the function gR{u)) whose behaviour can be understood by examining the 
equation (p.6|). 

Before we go ahead, we first discuss the asymptotic properties of the solution of 
the differential equation ( |3.6| ). Near the horizon the solutions are wave like and is 
given by 

(j)(r) {r — r j^) , r — )■ r+, (3.12) 
where Fh is given by the expansion 

F= (r-r+)F^ + ---, ^'^ = 4^- (3-13) 

Now at the boundary the two independent solutions are given by 

^ ^'^2(^) ~ r-^ r^oo, (3.14) 

(p[r) — — 7T—K2{t ) ~ constant, r — > oo. 

We will also need the fact that the differential equation given in ( p.6p can be obtained 
as the equations of motion of the following action 

= j^dvFr' (^|0'(r)p - ■ (3.15) 

3.2 Green's function in the w-plane 

In this subsection we will discuss the analytic properties of the function gR^u) defined 
in ( p.ll| ) in the complex cu-plane. 



No poles for Imoo > 0. 



Poles or divergences gR^uj) correspond to quasi-normal modes of the equation 
( p.6| ) |r3| . Quasi- normal modes are solutions to the differential equation ( |3.6| ) with 
the following boundary conditions. 

(f){r) {r — r^) , r — )• r+, (3.16) 
0(r) ~ r~'^, r — > oo. 

We will now prove that such quasi- normal modes do not exist. A simple intuitive 
reason that prohibits such quasi-normal modes is that, if quasi-normal modes with 



-7- 



Ima; > exists, they will correspond to instabilities due to the time dependence 
exp{—iut) in the mode given in (|3.5| ). We will now demonstrate such quasi-normal 
modes do not exist directly from the equation ( |3.6| ). All the coefficients in the 
equation ( p.6|) are real. Therefore if 0(r) is a quasi-normal mode with say complex 
frequency uj, then 0(r)* the complex conjugate is also a quasi-normal mode with 
frequency u*. This is seen by just taking the complex conjugate of the equation 
(^.61). We can now consider S^f, — S^jj = which by integration by parts and using 
equations of motion we obtain 

/•oo 3 

= Fr%<p*'<p-<p*(f)')\^^ + {u*^ -u^) I dr-\(P\\ (3.17) 



Here we have used the equations of motion of 0* in the first and the equation of 
motion of (f) in the second S^. Now from the conditions for the quasi-normal modes 
given in ( p.l6| ) and the fact that F ~ as r — >■ oo and F ~ (r — r/j)F/j as r — )■ r/^ we 
see that the first term in the above equation vanishes for Im(a;) > 0. Then we have 
the case that 

co^ = u;*^ (3.18) 

since the second term is positive definite. Thus we see that u has to be either 
purely real or purely imaginary, now with the condition Im(a;) > we find that u is 
restricted to be on the upper imaginary axis. 

Let us now assume that there exists a quasi normal modes on the upper imaginary 
axis with < 0. Then S^f, is positive definite. Substituting this solution in and 
integrating by parts we obtain 

S^ = Fr'<p*{r)<P'{r)\^^. (3.19) 



This vanishes from the behaviour in ( p.l6|) and the fact that Imoj > 0. But since 



is positive definite we must have = 0. Thus no quasi- normal modes exist in the 
upper half plane which implies no poles or divergences exist in the Green's function 
in this domain. 

No poles for cu real and to ^ 

For real oj, the following quantity 

W = (f)'(f)* - (f)'*(j), (3.20) 



is the Wronskian. From the differential equation in ( ^.6| ) we see that the Wronskian 
satisfies the equation 

W'+(^ + l)w = 0. (3.21) 
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The solution of this equation is given by 

W = ^. (3.22) 

We can determine the constant C since we know that if there exists a quasi-normal 
mode, the behaviour of the mode near the horizon is given by the ingoing boundary 



conditions given in ( |3.16| ). This fixes the C which determines the Wronskian to be 
given by 

W = (3.23) 

From here we see that the quantity r^FW does not vanish any where. We can 
now examine its behaviour as r — )■ oo for the quasi-normal mode. The boundary 
conditions for the quasi-normal mode given in ( ^.161 ) implies that — )■ as r — )■ oo. 
Using this fact in the definition given in ( |3.2(J| ) we conclude the quantity r^FW must 



vanish as r — oo which contradicts our previous conclusion. Thus quasi-normal 
modes and hence poles or divergences of the Greens's function do not exist on the 
real line for w 7^ 0. 

No poles for a; = 

We can determine the behaviour of the function as uj ^ Q explicitly from the 
differential equation (|3.6|) as follows. Let us define 

0'(^) 

m = (3.24) 

uj(j)[r) 



We can use the differential equation ( |3.6|) to obtain the equation satisfied by g. This 
is given by 

~g\r) + uj~g\r) ^ (L + l\ ^(^) + ^ = 0. (3.25) 



F r F"^ 



The ingoing boundary conditions at the horizon given in ( |3.16|) determines the bound- 
ary condition for g, this is given by 

~9 ~ ^ ~' V r ^ r+. (3.26) 

Fh{r - r+) 

The equation ( ^.251 ) can be easily solved in the a; — )■ limit. The solution which 
agrees with the boundary condition in ( p. 26 ) is given by 



g = -t^, cj^O. (3.27) 
This implies 

gR = oj lim r^Fg = -iur^, w 0. (3.28) 
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Thus we conclude that the Green's function has no poles as w — )■ 0. In fact as can 
be seen from the equation ( p.25|) , it admits an analytic power series expansion in uj 



at the origin. The leading behaviour is given above in ( |3.28| ). Using ( |3.25| ) one can 
easily set up a perturbative expansion in u and solve for gR^u) order by order in cu 
which gives rise to the power series expansion around the origin. 

No zero's for Imo; > 0, w 7^ 

Just as quasi-normal modes correspond to Dirichlet boundary conditions for the 
field at the boundary which leads to the boundary conditions given in ( p.l6| ) zeros 
of the Green's function correspond to Neumann boundary conditions for the field (p. 
This can be seen from the definition which implies that zeros of the Green's 

function correspond to the condition 

0(r) ~ (r-r+)"'^, r ^ r+, (3.29) 
r^F(j)'{r) — )■ 0, r 00. 

One might suspect that there might be singularities in which also gives zeros of the 



Green's function. However as discussed in [|13], a theorem of [0 guarantees that the 
solution is smooth with respect to a parameter u if the differential equation and the 
boundary conditions are smooth with respect to it. This ensures that there cannot 
be singularities in in the u plane. Now with the Neumann boundary conditions 
given in ( ^.29| ) we can repeat the steps of the discussion for the absence of poles or 
singularities in the Green's function. All the steps goes through since, here again the 
boundary conditions in (|3.29|) ensures that 

Fr^(0*'0 - 0>') = 0, r^r+,r^ 00. (3.30) 

Similarly for Imo; = 0, a; 7^ we have the limit 

WFr^ ^0, r ^ 00, (3.31) 

where W is the Wronskian. These conditions allows us to conclude that the Green's 
function does not vanish in the upper half u plane including the real axis, with the 
exception of a; = where we have shown that it vanishes and it admits an analytic 
power series expansion. 

Absence of branch cuts for Imu > 0^ 

Recall that the important properties of the Green's function is essentially con- 
tained in the function gR{u)) as defined in ( 3.11| ). In this equation 0(r) is the solution 



^ We would like to thank Aslioke Sen and the anonymous referee of this paper for useful comments 
on this section. 
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for the minimally coupled scalar equation in the background of (3^) with ingoing 



boundary ([4.14|) at the horizon. From (|3.11| ) we see that that we need to examine 



the behaviour of solution 0(r) and its derivative at the boundary. Note that the 
differential equation for the minimally coupled scalar and the boundary condition 
depends on the parameter uj and both these dependences are smooth with respect to 
this parameter u. Now the general theorem in ^ ensures that the local Forbenius 
expansion of the solution is smooth with respect to the parameters of the differential 
equation if both the equation and the boundary condition are smooth with respect 
to the parameter. Hence the Forbenius expansion of (f) and its radial derivative at 
the boundary are smooth with respect to the parameter uj. Now using this theorem 
it is easy to see that the n*'^ order derivative of the gn^u) with respect to u is smooth 
in the upper half u plane. Since and its derivative are smooth with respect to uj 
we see from ( p.lip the only possible locations of singularities of the n^^ order deriva- 
tive of the gniuj) are the zeros of 0, that is where it satisfies the quasi-normal mode 
boundary conditions given in ( ^3.161 ). However we have shown earlier that there are 
no-quasi-normal modes in the upper half-plane, this assures that arbitrary deriva- 
tives of gniuj) exists in the upper half plane. This implies that the Green's function 
does not have any branch cuts in the upper half-plane. 

Thus we have concluded that the Green's function satisfies property 1 for the 
derivation of the sum rule. We now need to ensure that property 2 is satisfied. 

Behaviour as a; — t- oo 

To obtain the behaviour of the Green's function at large u we first rewrite the 
differential equation given in (|3.6| ) by defining a dimensionless variables 

z = — , i\ = — u. (3.32) 

r r_|_ 

The equation reduces to 

- + z')<P'{z) - = 0' (3-33) 

where f{z) = l — z^. Note that for convenience we have gone over to the Euclidean 
frequency labelled by A. Let us re-state the boundary conditions in terms of these 
coordinates. The ingoing boundary condition reduces to 

A 



l-z)4, z^l (3.34) 



^In fact this theorem has been used in the identification of the location of the poles of the Green's 



function with the quasi-normal modes of in 1 13 . Since in general the poles of the Green's function 



can also arise from the singularities of 4>' . However this theorem guarantees that </>' is smooth with 
respect to w hence poles can arise only from the zeros of </> which correspond to its quasi-normal 
modes. 
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and the Green's function is given by 



4 



9R{X) = —^]im-!-. (3.35) 

Our aim now is to find the behaviour of this function as A — )■ oo. For this purpose 
we first rescale the co-ordinates as 

y = Xz. (3.36) 
The differential equation given in ( p.33| ) can be written as 



- + J,)ny) - -j^.m = 0, (3.37) 



4 

where f = 1 — and derivatives are with respect to y. We can solve the equation 
using the Forbenius expansion at around y = X. This results in following series 



oo 



The coefficients aj can be determined by recursion relations on substituting the 
Forbenius expansion in the differential equation ( p.37| ). Now the differential equation 
in ( p.37|) also admits an expansion at around y ^ 0. This expansion can be organized 



as a systematic expansion in powers of j as follows: We first expand ( |3.37| ) in powers 
of 1/A, this results in 

<P"{y) - + 4|J + + ■ ■ .)<p'{y) - (1 + 2|J + 3|^ + ■ ■ ^^(y) = 0. (3.39) 

One important fact to point out in this expansion is that at A strictly infinity, the 
equation reduces to that of the minimally coupled equation in the background of 
pure AdS^, that is the background corresponding to zero temperature. To solve this 
equation perturbatively in ^ we define the quantity 

(3.40, 

From the expansion of the equation in ( |3.39| ), we see that the we can expand g as 

9 = 9o{y) + ^9i{y) + ^92{y) + ■ ■ ■ (3.41) 



CXD ^ 

n=0 
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We can obtain the equations satisfied by each gn by substituting the expansion above 
and matching orders in The first few equations are given by 

3 

9o{y) + gliy) - -9o{y) -1 = 0, 

g[{y) + (2go{y) - gi{y) - Ay'go{y) - 2y^ = 0, (3.42) 
M + (^goiv) - ^) g2iy) + gl{y) - ^y^gi{y) - ^y'g^{y) - = 0. 

The advantage of casting the equations in this form is that the apart form the first 
equation, the subsequent equations are first order equations. We will now illustrate 
how to obtain the two independent solutions. There are two solutions to the first 
equation, these are given by 

Let us first find the expansion around the first solution, for that we have to solve 
for gi, this is given by 



^2/^^(1-^)^-^. (3.44) 

The term proportional to Ci is the homogeneous solution. The constant Ci is set 
to zero by requiring that the asymptotics of the solution set by g^ be unchanged in 
the presence of gi. Note that if Ci ^ 0, gi grows exponentially when y ^ oo. This 
condition is imposed on all equations. Thus we have 



,(i)..A_o..3 , y\, ^ 



^r(y) = 2y^ + ^(l-^), (3.45) 

"2 



= -^2/'-^|/' + 0(/,/log(y)). 
Similarly the solution of g2 is given by 

g2{y) = rc^y[-(^«)^(y)+4y=^^«(y)+4y^5?^(y) + 3y«](yK,^(y)). (3.46) 

yKiiy) Jo 

Prom the expansion of the functions involved, it can be shown that 512^^ admits the 
following expansion near the origin 

g?{y) - + y' + (log(y) + log(2) + ^)y'' + ■■■ . (3.47) 
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An simple examination of the equation for Qn shows that the leading term in the 
expansion of gn^ around the origin is given by 

m>7, forn>2. 



Now that we have Qn^ the first solution (p^^\y)) is given by 
0«(y) = exp(/ dyJ2{^), 



n=0 



(3.48) 



(3.49) 



y'My) 1 



1 1 
^A^ 2 



Jo 



dy9i\y) 



+ 



Note that (p^^^ has the following behaviour as ?/ — 



constant. 



(3.50) 



From (|3.45|) we see that the ^ coefficient in the expansion given in ( p.49| ) behaves as 
y'^ near the origin. The equation ( p.48|) shows that the remaining terms are further 
suppressed as y — )■ 0. 

A similar construction can be done starting with the seed gl^\ The first ^ 
expansion around the second solution is given by 

ry 



9?\y) 



1 



y^iy) Jo 



dy[4y'gf\y) + 2y\yl',{y)), 



^f)(y) = 2y^ + ^(l-l|). 



5 



(3.51) 



Note that we have set the constant corresponding to the homogeneous solution to 

(2) 

zero. Similarly the solution of g2 is given by 



9i'\y) 



1 



dy[-{g'i'^)\y) + Ay^g^'iy) + ^y'9^''{y) + 3i/«](y/|(y)). (3.52) 



,3^(2) 



7^(2), 



yiKy) Jo 

From the expansions of the functions involved it can be shown that 

9?\y)-y'', y^o. 

Similarly in general it can be shown that 

9l^\y)-y'^, m>7, forn>2, y -> 0. 
Thus we obtain the second solution 



^('){y)=y^I,{y)ll + - I dyg^'iy) 



.(2)/ 



+ 



1 1 

A8 2 



Jo 

dyg?\y)f+ I" 9?\y) 



(3.53) 



(3.54) 



(3.55) 



+ 
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It will be useful to note that 0*^^^ has the following behaviour as ?/ — )■ 0. 



0. (3.56) 

Note that from (|3.51|) and (|3.54|) we see that the coefficients of the ^ expansion are 
further suppressed as ?/ — )■ 0. 

Once the solutions 0*^^^ have been constructed, then the solution which 
matches the ingoing boundary condition at the horizon is obtained as follows: Con- 
sidering the linear combination 

(/.(?/) = A(A)0«(l/) + i?(A)0(2)(y). (3.57) 

The coefficients A, B which depend on A are obtained by matching the solution (j){y) 
in (|3.57|) and that obtained by the Forbenious expansion aX y = \ given in (|3.38|) . 



The expansion in ( p.38|) is asymptotically expanded to ?/ — )■ 0, while the ^ expansion 
in (|3.49|) is expanded to ?/ — ?■ oo and the function and the derivatives are matched to 
determine A and B. Though we do not explicitly require A and B as functions of A, 
we can argue that 

fi(A)^0, A(A) ^ ^, y^oo. (3.58) 



The reason for this is that when A strictly oo, the equation ( p. 391) reduces to that of 



the zero temperature case as we have observed earlier. In this situation the solution 
which is finite at y — is given by ^y'^K2{y) ^. Thus we arrive at ( 3.58 ). 



Now let us study the implications of the properties of the solutions of the differ- 
ential equation (|3.37|) have on the retarded Green's function Qr^uj). In terms of the 
scaled variables it is given by 

a (X) - hm -4(A)0W'(y) + i^(A)0(^)'(y) 



Now from the behaviour of the solutions as ?/ — )■ given in ( |3.45[ ), ( |3.48| ), ( |3.51| ), 



( p.54|) and the discussions below (|3.49|) , ( p.55|) we see that the y ^ limit in ( p.59| ) 



reduces to 



L2 yy^o y3 ^0 vy; ^ ^ ^ ^^^^ j 

In the second line of the above equations we have substituted the expansion of g^^ (y) 
near the origin which is given in ( p.45|) . An important point to note is that the above 



^The 5 is a normalization such that the solution ^y^K2{y) ^ 1 as y — > 0. 
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result is true for all values of A or frequency. We can further take the limit A — )■ oo 
and using (|3.58|) we are left with 

hm gniX) = flim^A^^« - |) . (3.61) 

As we have observed the leading contribution g^"^ is identical to the zero temperature 
case, it is divergent as A — )■ oo. From the expression of g^^ in ( |3.43| ) we see that 
this term also has l/y^ and log(?/) divergences as ?/ — )■ and needs to be regulated 
with appropriate counter terms. The next order correction in this limit is the finite 
contribution from the constant piece involving gi. Thus the behaviour of the Green's 
function does not satisfy property 2. 

To regulate the Green's function so that it satisfies property 2 we follow |Q and 
consider 

5Gr{uj) = Gn{u, T) - Gr{u, 0) + + P- (3-62) 

Let us recall the definition of each of the terms in the above expression. 

Gr{uJ, T) = - gR{u) - P + G'countcr(w, T). (3.63) 

Here we have just substituted the expressions for Gr{u) and Gcontact in the expression 
for the Green's function given in (|3.^). The T = Green's function is given by 



J\t2 »,4 -I 

Gr{cO, 0) = ^YT^J^ -^'9i'\y) + Counter (CU), (3.64) 
- -OJ hni h Gcounter(w,0). 



Note that there is no contact term in this case since at T = 0, the pressure vanishes. 
We now point out that we have the following equality 

Gconntcri^,^) = G counter ^T) . (3.65) 

The reason is clear because the and log(r) divergences occur only in the term g^"^ 
and this term is common both for T 7^ and T = 0. Thus the counter term to cancel 
these divergences must be identical. We also can infer that they are proportional to 
o;^. Now let us examine the behaviour of 5Gr{uj) as w — )■ 00. Substituting all the 
equations we see that 

lim 5Gr{uj) = lim ^ ^ 0. (3.66) 

^^00 ' A^oo Svr^L^ 4L2 v4(A) 

Thus bGRibj) satisfies the property 2. What has been essentially done by the con- 
struction in (|3.62|) is the divergent and constant pieces of the Green's function as 
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u ^ CO are subtracted out. Now that 6Gr{u!) satisfies property 2. we can apply 
Cauchy's theorem and obtain the sum rule. For later purposes note that 



Im 6Gr{co) = Im Gr{uj, T) - Im Gr{{uj, 0), (3.67) 

= p(^^,^) -p(^,o), 



since the remaining expressions in (|3.62|) are real. 



3.3 The sum rule 

Now that property 1 as well as property 2 are satisfied by the function 5Gr{uj), we 
can evaluate the LHS of the sum rule. This is given by 

5Gr{S)) = Gr{0, T) - G^(0, 0) + ^^7^3 + P. (3.68) 

Now we also have 

Gk(0, T) = -^-^^^(0) - P, (3.69) 

= ~P + G'counter(0, T), 

Gh(0,0) =0. 

Here we have used that fact that Qr^Q) = and the fact that Gcounter vanishes at 
a; = since as we have seen earlier that they are proportional to w^. Using the 
equations ( |3.69|) in ( p.68|) we obtain 



7V^ r1 6 2 

^-Gh(0) = = (3.70) 



Now using (|3.67|) and the fact that 6Gr{u) satisfies both property 1 and property 2 



we obtain the shear sum rule derived in 

-e = - / — {p{uj) - Pt=o(w)) • (3.71) 

Note that the LHS of the sum rule essentially originates from the constant term in 
the function gR^u) as u ^ oo. 

3.4 Sum rule from OPE 

From the explicit derivation of the sum rule, we see that the LHS side essentially 
depends on the high frequency behaviour of the Green's function. In fact it is just 
proportional to the constant term in the high energy behaviour of the Green's func- 
tion. Coefficients in the Operator product expansion of the stress tensor contains the 
information of the short distance or high energy behaviour of the Green's function. 
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Thus we should be able to obtain the sum rule from the OPE of the stress tensor . 
This OPE is given by 



T^u{x)TpM ~ C-T^ + Ki^paap{x)T^p{Q) + ■ ■ ■ . (3.72) 



A contains various Lorentz structures which can be found in |T6|. The important 
property of them which is seen form conformal invariance is that they all scales like 
Now we take the Fourier transform of the above OPE and set g = and 
w — 7- oo. By a simple scaling analysis it can be seen that the term proportional 
to ^ scales like cu^. This divergence is the same divergence seen in the gravity 
calculation due to the terms which are identical to the zero temperature case. The 
term proportional to A gives rise to a constant independent of u. Thus the constant 
contribution to the Green's function as w — )■ oo arises due to the one point function 
of the stress tensor in the thermal ensemble. Now evaluating 

lim {Gn{uj, T) - Gr{uj, 0)) = -P - \t. (3.73) 

We have used the result obtained in ||^ but in the Minkowski space and the convention 
of From hydrodynamics since the Green's function is defined as the response of 



the changes in the stress tensor on perturbing the metric, the zero frequency value 
of the Green's function is just given by negative of the pressure. 



Gr{Q,T) = -P, Gr{lo,Q)=Q. (3.74) 



On defining 



5Gn{ujT) = Gr{uj, T) - Gn{uj, 0) + P + (3.75) 

5 

we see that the LHS side of the sum rule is just the difference between the con- 
stant terms in the high frequency behaviour and the zero frequency behaviour of the 
Green's function. This results in the sum rule given in ( p.71|) . 



4. The shear sum rule for the R-charged D3-brane 

In this section we examine the gravity dual of A/" = 4 Yang-Mills at finite chemi- 
cal potential and finite temperature and re-derive the sum rule. Since there are 3 
R-charges corresponding to the Cartan's of S'0(6), it is possible to turn on 3 indepen- 
dent chemical potentials. The gravity dual of this system is given by the R-charged 



black hole of Behrndt, Cvetic and Sabra [17|. Using the differential equation of the 



massless minimally coupled scalar in this background, we obtain the retarded Green's 
function of the T^y component of the stress tensor. Examining this differential equa- 
tion and the method discussed for the uncharged D3-brane, we will show that the 
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regulated Green's function satisfies both property 1 and property 2 which is neces- 
sary for the derivation of the sum rule. We show that the sum rule is corrected by 
terms which depend on the 3 chemical potentials. We explain these additional terms 
in the LHS of the sum rule due to the presence of scalars in the OPE of the stress 
tensor. The procedure to demonstrate property 1 and property 2 of the retarded 
Green's function parallels that of the uncharged D3-brane, therefore we will be brief 
and highhght only the differences. 

4.1 Green's function from gravity 

We begin with the metric for the R-charged D3-brane with all the three charges 
turned on which is given by 

dsl = (Z^W / dt' + 7^V3 + dy' + dz') + n'/^^du' 

U U ^ ' 4:JU^ 

3 

f(u)^n(u)-u^ll{l + h), Hi^l + hu, h=^, To^^, 

i=l H 

u^^, = (i + M(i + M(i + M- (4-1) 

The scalar fields and the gauge fields in this background are given by 



^^ = ^' ^^-T^V h-V^.Lflil + k.r. (4.2) 

i=l 

The above metric is a solution of the equations of motion of the STU-model given 
by action 

S-^J d'xV^C, (4.3) 



where F'^^^^i = 1,2,3 are the field-strengths for the three Abelian gauge fields. The 
three scalar fields X*'s are subject to the constraint X^X'^X^ — 1. The metric on 
the scalar manifold is given by 

G,. = idiag{(X^)-^(X^)-^(X^)-}. (4.4) 

The scalar potential is given by 

V = 2(i, + i, + 3L). (4.5) 
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It is useful to recall the following thermodynamic data of this black hole |18|. The 
Hawking temperature Th, entropy density s , energy density e, pressure P, charge 
densities pi and the conjugate chemical potentials pi are given by 

2 + ki + k2 + h- kikih -k'^N'^T^ , , , 

^^ = o u n , .7°' ^ = ^^11(1 + ^^^)''^ (4-6) 



2v/(rTW+W+^ 2 



— ^ + ^^)' ^ = 11(1 + 

i=l i=l 
3 m /TTT- 3 



The thermodynamical stability condition of this black hole is given by 

2-ki-k2-k3 + kik2ks > 0. (4.7) 

As before, to study the retarded correlator Gr{u, T) we must examine the equa- 
tion of motion of a massless minimal coupled scalar in this background. We consider 
the perturbation 

5g,y = <p{r)'H'/'^e-'^\ (4.8) 
whose equations of motion in radial co-ordinate r is given by [|18 



5.^0+ + ^0 = 0, (4.9) 



where again F is given by 



and / is defined in (f4.1|) . This equation is similar to the uncharged case except for 
the presence of H as coefficient of cu^ and the fact that / depends on the constants 
fcj. To obtain the retarded Green's function, we must impose in-going boundary 
conditions at the horizon r+ and obtain at the boundary r — )■ cxd. Then the 
Green's function is given by the expression in (|3.^ ). Gcontact is the contribution of 
the contact terms obtained from the on shell effective action. For the charged case, 
this can be evaluated by using following action 

Here the second term is the Gibbons- Hawking boundary term while the third term 
is the term required to make the action finite in the limit r — )■ oo [11^] . It is given by 

-1/1/3 3 
i=l 
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Using this, it can be shown that the contact terms are again obtained from the on 
shell action given in (|3.9|) [ p!8[| ^. Therefore we obtain 



Ccontact — —P- (4-13) 

There is an alternative indirect method to infer that the contact term is given by 
—P. Contact terms do not depend on frequency, therefore one can obtain them by 
taking the zero frequency limit. It will be shown that G{u,T) — )■ in the w — )■ 
limit, further more the Gcountcr as in the case of the uncharged case behave as u^. 
Thus in the zero frequency limit the only contribution to the Green's function is from 
Gcontact- Now from hydrodynamics, the Greens' function is the response of the stress 
tensor T^y to the metric fluctuation 6gxy. At zero frequency, this is just —P. Thus 
we have Gcontact = 

Before we begin our analysis, we will need the behaviour of the solutions near 
the horizon and the boundary. The behaviour of the solutions near the horizon is 
given by 

0(r) ~ (r - r+)^^"'^, (4.14) 

where a can be written as 

« = ^Vi^ + ki)il + k2)il + h) (4.15) 

F = {r- rh)Fh + • • • , = ^{2 + + k2 + h - kik^k^). 

As r — )■ oo, the metric asymptotes to AdS^, therefore the behaviour of (f){r) is the 
same as in the uncharged case and is given by ( p.l4|) . Finally the equations for 0(r) 
can be obtained by variation of the following action 

S, = fdrFr' (^|0'(r)|2 - ^|0(r)p) . (4.16) 

The only change in this action compared to the uncharged case in ( |3.15|) is the 
presence of "H in the term proportional to the frequency. 

4.2 Green's function in the tj-plane 

In this subsection we will discuss the analytic properties of the function 

Fr^d)' 

gniu) = lim — (4.17) 

r— >oo (7) 

in the complex u plane. The discussion closely follows that of the uncharged case. 
The difference arises mainly in the behaviour of the w — )■ oo limit which we will 
highlight. 



*See equation 4.24 of 
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No poles for Imo; > 



Here the discussion is identical to that of the uncharged case. Poles or divergences 
in gn^oj) correspond to quasi-normal modes of the equation (|3.6| ). Quasi-normal 
modes are solutions to the equation ( p.6| ) with the following boundary conditions 



(r) ~ (r - r+)e-*"'^, r ^ r+, (4.18) 



fr) ~ r ^, r — > oo. 



One can follow the same arguments as in the case of the uncharged situation to show 
that such modes do not exist for Imu > 0. We need to employ the action given in 
( [4. 161 ) for this purpose. The fact that Ti > in the domain r+ < r < oo also ensures 



that the action evaluated on the quasi-normal mode is positive definite. This 
enables us to repeat the same arguments as in the case of zero chemical potential. 

No poles of oj real and oj ^ 

Here again the same reasoning developed for the uncharged case earlier using 
the Wronskian, ensures that quasi-normal modes do not exist on the real line u ^ 0. 
This implies that there are no poles or divergences in gR{uj) in this domain. 

No poles for a; = 

The behaviour of the function ^f^j as w — )■ can be determined from the differ- 
ential equation ( [4.9|) as follows. We define 

9{r) = ^r-y (4.19) 

uj(p[r) 

Using the equation ( [4. 9] ) we find the g{r) satisfies the equation 

9+-f+[^ + l)~9ir) + ^ = 0. (4.20) 

The ingoing boundary condition implies that g satisfies 

_,^(l + fc^)(l + A:2)(l + fc3) ^ 
Fh{r - r+) 

The solution of ([4.20|) in the u —> limit which satisfies the above boundary condition 
is given by 

g = -i-^ ^, w ^ 0. 4.22 

This implies that 



lim gR{Lo) = cu lim r^Fg = -icorlJil + fci)(l + k2){l + h). (4.23) 

UJ—^0 r—>-oo 
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Thus gRioj) is proportional to w as w — )■ 0. Examining the equation (|4.20|) , it can be 
seen that it admits a power series expansion in u around the origin. 



No zeros for Imo; > 0, a; 7^ 

As argued for the uncharged case, zeros or points at which the Green's function 
vanish correspond to the boundary condition 

0(r) ~ (r - r+)-^°'^, r ^ r+, (4.24) 
r^F(j)'{r) — )■ , r — )■ 00. 

Using the same reasoning discussed for the uncharged situation, it is easy to see that 
there are no modes which have the above boundary condition for Imw > 0, 7^ 0. 
This imphes that there are no zeros or points at which the Green's function vanishes 
in this domain. 

No branch cuts for Imw > 

From the discussion above, we have seen that the Green's function does not have 
any singularities, nor does it vanish any where in the upper half plane except at the 
origin. At the origin it admits an analytic power series expansion. As we have argued 



for the uncharged D3-brane case the general theorem in (||T^) ensures that the local 
Forbenius expansion of the solution is smooth with respect to the parameter u since 
the differential equation and the boundary condition are both smooth with respect 



to this parameter u. Now from ( [4.17| ) it is easy to see that the rfi^ order derivatives 



of the Green's function with respect to the parameter u exist. This is because the 
singularities for the n^^ order derivative can arise only from the zeros of which 
are its quasi-normal modes. We have already shown that there are no quasi-normal 
modes in the upper half-plane, hence we conclude that arbitrary derivatives of the 
Green's function with respect to u exists in the upper half plane. Thus the Green's 
function does not have any branch cuts in the upper half-plane 

Therefore the expansion at w = can be extended to the entire upper half plane 
and the function is analytic in the upper half plane. Thus the Green's function 
satisfies the property 1 for the derivation of the sum rule. 

Behaviour as a; — )■ 00 

To obtain the asymptotic behaviour of the Green's function it is again convenient 
to work in the scaled variables 

'''' i\ = —u. (4.25) 



+ 
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We then expand the equation (|4.9|) in powers of j. For this it is convenient to express 
the background metric given in ( [4.1|) in Fefferman-Graham coordinates. The metric 
for the charged D3-brane in Fefferman-Graham coordinates is given by 



ds^ = ^{ (-1 + J2ai{zy)dt^ + (1 + J2bi{z^y)idx^) + dz^\. (4.26) 

^ \ 2=2 1=2 / 

The details of the transformation to these co-ordinates are given in the appendix. 
Note that in this expansion the term ai, bi is missing. Here the boundary is at 2; = 0. 
It will be seen that the coefficients a2, 62 are crucial to determine the asymptotic 
behaviour of the Green's function. These are evaluated in ( |A.5| ). The equation for 
the minimally coupled scalar in this co-ordinates is given by 



\ 1=2 J \ i=2 



z'y U = 0, (4.27) 



where aj,6j are related to the coefficients aj,6j. Again the leading coefficients are 
important for our later use which are related by 

02 = 2(02-862), 62 = 02. (4.28) 

We now rescale the co-ordinate z by defining 

y = \z. (4.29) 

This leads to the equation for the minimally coupled scalar as a ^ expansion. 

T / 00 2 \ CO 2 

- ^ 3 + E^^(fi)^ ^'(y) - (1 + T.~^^i^y^<Piy) = 0- (4-30) 

y \ i=2 J i=2 

Now this expansion is similar to that obtained for the un-charged situation given in 
( ^.491 ). We can therefore use the same method to solve it order by order in 1/A. We 
define 

From the equation given in ( [4.30|) we see that the g admits an expansion of the form 

11 , , 

9 = 90 + -^91 + j^92 + ■■■■ (4.32) 

The equations satisfied by these can be obtain from (^4.30| ) and the first few orders 
are given by 

^7o + ^?o-^-l = 0, (4.33) 

g[ + (^2go -fjai- a2y^9o - kv' = 0, (4.34) 
3 

y 



g'^+[2go--]g2- a^y'go - hy' = 0. (4.35) 
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We will now follow the same procedure as done for the uncharged case to obtain two 
independent solutions to ( [4.30|) . The zeroth order solution is same as the uncharged 
case since in the A — > oo limit, the equation ( |4.30| ) reduces to the A — )■ cxd hmit 
of ( |3.49|) . This is expected since asymptotically the metric ( [4.1|) reduces to AdS^. 
Therefore we have 

Let us first find the j expansion around the solution gl^\ The first order correction 
is given by 

We set ci = 0, so as not to change the asymptotics of the solution at ?/ — )■ oo. On 
performing the integral we obtain 



10^ ^ 



^7r^(z/) = 7fr + 7^1/^(1-^), (4.38) 



2 ^b2jy' -—b2y' + 0{y\y'\og{y)). 

Similarly, the solution to g2 is given by 

9?{y) = ^jj^^ I {asv'go + hy') yKl{y). (4.39) 

From the expansions of the functions involved, one sees that 

g'i\y) = y'' + 0{y\yHogy). (4.40) 

Similarly from the structure of the functions in the j expansion, it is easy to see that 

gn{y) -2/", rn> 7, for n > 2, y -> 0. (4.41) 

Using this procedure it is possible to construct all the gn^'?>. Form this we find that 
the first solution of the differential equation is given by 

0(1) (y) = eMjyy (^0 + + f + ■ ■ ■) , (4.42) 

= y^K2{y) [l + Y^ dygiiy) + ^ ^ dyg2{y) + ■ ■ . 

Note that 

0(1) ~ constant, y ^ 0. (4.43) 
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From ( [4.3S ) we see that ^ term in (|]4|) goes as 7/^ near the origin. The higher 
order terms are further suppressed as y — )• 0. A similar construction can be done 
starting with the seed . The first order term about this seed is 

9f\y) = dy (a^y^go + h^y') , (4.44) 



9?\y) = ^jjf^^ I {dsv'go + hy') yll{y). (4.45) 



0,2-1. &2 c: / 

Similarly g^\y) can be written as 

ry 

yiliy) 

Examining this solution, we see that the behaviour of this function as ?/ — )■ is given 
by 

(2/)~y' + 0(/), y^O. (4.46) 
From the general form of the equations for gn-, it can be seen that 

9n\y)--y'^^ m>5 for n>2, y ^ 0. (4.47) 
Thus one obtains the solution 

'P^'Ky) = y'hiy) (i + Y,J\i{y) + [ a^iy) + ■ (4.48) 

It is important to observe that 

(t>^^\y)--y\ 2/->0, (4.49) 

which is same as in the case of the uncharged case. Again from ( [4.44| ) and ( |4.47| ), 
the terms in the j expansion are suppressed as y — )■ 0. 

Once the two independent solutions have been constructed in this manner, we 
can write the solutions which obey the boundary conditions as 

0(y) = A(A)0«(y) + i?(A)0(2). (4.50) 

The coefficients A, B depend on A and can be obtained by matching 0(?/) and 
its derivative with the near solution. The near solution is obtained in terms of a 
Forbenius power series expansion at the horizon. The solution is of the form 

0=(l-f)'^(^E«,(l-fyj- (4.51) 

Now though we do not need the detail dependence of the coefficient B on A, it is 
clear that using the same argument as in the uncharged situation we find 

B{\) 0, A{X) -> ^, as A oo. (4.52) 
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The reason is that we know that at A — )■ oo, the equation reduces to that of the 
T = hmit. In this situation the solution which is finite at y — )■ oo is ^y'^ K2{y) . 

With this input and using the same steps followed in the uncharged situation we 
find the Green's function to be 



In this we have substituted for 02, &2 from (f4.28|) . Again, the above result is valid for 



all values of A. Taking the A — )■ 00 limit we are left with 

lim 9RiX) = I lim ^^^^^M -ha, + 5b,) \ , (4.54) 



where we have used (^4.521). As we have seen earlier the A^ contribution is identical 



to the zero temperature case. There is also a constant term due to the first order 
correction g^^\ Thus the Green's function does not satisfy property 2 and needs to 
be regulated. 

As before, to regulate the Green's function we consider 

5Gr{u) = Gniuj, T) - Gr{uj, 0) + f ^ («2 + 562) 1 + P- (4.55) 

Now using all our previous results and following the same steps as in the case of the 
uncharged situation we have 

6Gr{uj) ^0, a; ^ 00. (4.56) 

We have essentially subtracted the divergent and the constant term in Gr{uj) so that 
property 2 is true on 5Gr{uj) which can then be used to obtain the sum rule. Note 
that for the charged D3-brane case the equation ( |3.67|) still holds since all the terms 
subtracted to construct 6Gr{uj) are still real. 

4.3 The sum rule 

We can now use Cauchy's theorem on the function 6Gr{uj) and derive the sum rule. 
The sum rule for the charged D3-brane then is given by 

SGniO) = r ^ {p{u) - pT=o(u;)) , (4.57) 
J-00 ^ 

where 

TV^ r*^ /3 \ 
^^^(°) = 8^Z^(5^"^ + ^^^V- ^^-^^^ 
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We can now substitute the values of 02, 62 from ( A. 5 ). We obtain 



^t2 4 o 

'^^^^^^ = s^^Z^s ^^^^ + ^^^^^ + ^'^^^ + ^^-^^^ 

{(^1 - ^2)' + (^1 - ^3)' + {k2 - A:3)'}) 

Here we have rewritten the first term in terms of the energy density using its expres- 
sion given in ( |4.ti| ). Note that there is a modification of the sum rule which involves 
the chemical potential. In the next section we will show that this modification is due 
to the expectation values of the scalars in the charged D3-brane background. As a 
simple check of this hypothesis, note that for D3-brane with equal R-charge in all 
the 3 Cartan directions ki = k2 = k^, the additional term vanishes. 

4.4 Sum rule from OPE 

We have seen that the term on the LHS of the sum rule is essential due to the high 
frequency behaviour of the retarded Green's function. In the uncharged situation, 
this high frequency behaviour can be captured by the OPE of the T^y component 
of the stress tensor. Therefore if the sum rule has to be modified there must be 
additional terms in the OPE which contribute at finite chemical potential. Let us 
now re-examine the OPE of the stress tensor which is given by 

T^.(x)T,,(0) ~ cj-^ + i^.paa/3(x)T„^(0) + 5;,^,(x)a(0). (4.60) 

Here we have included additional terms in the OPE which arise if there are scalars 
denoted by Oa in the theory and the three point function {T^uTpaOa) is non-zero. 
From a simple scaling analysis it can be shown that operators of dimension lower 
than 4 in the OPE contribute to divergences as a; — ?■ 00. Operators of dimension 
greater than 4 are irrelevant and operators of dimension 4 contribute to the finite 
terms responsible for the sum rule. Consider an operator Oa of dimension A in OPE, 
the structure of this term will be of the form 

B;,^,ix)OM, (4.61) 

where a refers both to the operator as well as tensor indices if any. By conformal 
invariance we have the following scaling property 

B;,^^{Ax) = K^~'B;^^^{x). (4.62) 

Now, on taking the Fourier transform of this term we obtain 

j rf^xe— *i?^,^,(x)a(0) = u'-'' I d^£e-*i?^,^,(x)a(0). (4.63) 
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Here we have rescaled all coordinates by replacing x = ^. From ( [4.63 ) we see that 



the terms diverge for A < 4 and are irrelevant for A > 4. The finite terms in the 
OPE are therefore coming from operators of dimension A = 4. 

We have seen from the discussion of the uncharged D3-branes, the term propor- 
tional to the stress tensor contributes to term |e in the sum rule. Thus the remaining 
terms in ( |4.59| ) must be due to the expectation values of certain dimension 4 opera- 



tors. From the structure of the extra terms in ([4.59|) , we see that they vanish when 
all the charges are equal. At this special point one can see from the solution ( ^4.2| ), 
the scalars become trivial while the gauge field still remains non-zero. Thus it must 
be the case that operators dual to the scalars are responsible for the additional terms 
in the sum rule. To find the conformal dimensions of the operators dual to these 
scalars we evaluate their masses. Due to the constraint X1X2X3 = 1, there are two 
independent scalars. We parametrize Xj's in terms of these independent scalars by 

Substituting these redefinitions in the action ( [4.3| ) and expanding the scalar potential 

v-T±i. («5) 

i=l 

to quadratic order in 0j it can be shown that the masses of both 0i, 02 is given by 

m^L^ = -4. (4.66) 

From the mass- dimension relation 



A(A - 4) = m^L^, (4.67) 

we see that the fields 0j corresponds to operators of dimension 2 in the field theory 
and saturate the Breitenlohner-Freedman bound. In the field theory these operators 
correspond to two linear combinations of the following chiral primaries 

Tr(XX), Tr(FF), Tr(ZZ), (4.68) 

where X, Y, Z are 3 complex fields constructed out of the 6 scalars in A/" = 4 Yang- 
Mills. These are the three scalars which are un-charged under the Cartans of 5*0(6). 
Note that the combination Tt{XX) + Ty{YY) + Tt{ZZ) is the Konishi scalar and 
therefore not a chiral primary. Let us denote the two chiral primaries dual to the 
field 01,02 as Oi and O2 respectively. 
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From the gravity background it is easy to read out the expectation values of 
the operators Oi and O2 at finite chemical potential. Since these fields saturate the 
Breitenlohner-Freedman bound their behaviour at the boundary is given by 

0.(r)~L(a.^ + Al + ...). (4.69) 

Here we have introduced factors of L so as to ensure that aj, Pi have the appropriate 



dimensions. The boundary term for the scalars given in ([4.12|) ensures that ctj is set 



to zero at the boundary. Then the expectation value of the operator dual to (pi is 



obtained by substituting the redefinitions ( [4.64| ) and in the background value of Xi 
and reading out the value of g^^A pO[. See for earlier work. They are 

given by 

^^'^ - 8^^^ - 8^Z^ — 7! — ' ^^-^^^ 

^ ^' M^L^'^ 87r2L4 72 ■ 

Note that these expectation values have mass dimension 2 since the operators Oi, O2 
have mass dimension 2. They are also constant in space time. 

A natural candidate for the operators of dimension 4 which occur in the OPE 
and which can contribute to the sum rule are the squares of the operators ^ 

(^1 = ^0^, 0, = ^0l 0, = ^0,0,. (4.71) 

We have defined these operators with the factor of so that the two point function 
of all these operators are normalized as A^^. 

We first show that the OPE of the stress tensor can give rise to these operators. 
The operators Oi and O2 are chiral primaries, and they will contain say the operator 
Tr(ZZ). Thus the operator Oi and O2 will contain the operator 

V = ^mZZ)f (4.72) 

Therefore we evaluate the following three point function in the free field limit 

(T,,(x^)T,,(x2)P(x3)) = ^(Tr(9,iZ9,iZ)Tr(9,2Z9,2F)Tr(ZZ)Tr(ZZ))(4.73) 

To obtain this equality we have written out only the terms which contribute to the 
correlator in the free field limit. We have also ignored numerical proportionality 



^The operators Oi do not occur directly in the OPE since the 3 point function {T^^Tp^Oi) 
vanishes in gravity. 



- 30 - 



constants. Note that the two point function of the stress tensor is normahzed to 
unity. By applying Wick's contractions we obtain 

(T.,(.W)n^')} = /V - ^7^'; ^><-^ : - (4,74) 

00 tjC ijC ijC 

Therefore from the free field calculation it is evident that 

{T^yT^yV) ^ (4.75) 

Any operator O whose three point function (TTO) is non-zero will appear in the 
OPE of the stress tensor. In fact they will be more relevant than the stress tensor if 
the operators have dimension A < 4, and as relevant if the operator has dimension 
A = 4, see |T6| below equation (6.47). Thus even in the free field limit the operators 



considered can arise in the OPE of the stress tensor. 

Now we will show that the expectation values of the operators Oi and O2 is 
responsible for the modification of the sum rule. For this we evaluate the following 

{O^) + (O,) = + {0,f) (4.76) 

= ^ {{ki - hf + - hf + {k, - hf] . 

The first equality is due to large factorization. In the second line we have substi- 
tuted the expectation values form ( [4. 701) and re-expressed r_|_ in terms of Tq. Note 
that the final line of ( |4.76| ) is proportional to the additional term in the sum rule 



( [4.59| ) with the exception of the factor of A^. This fact is because the structure 
constant -B^^p for these operators scales as in the large A^ limit as can be seen 
from the normalization of our operators. This additional factor of together with 
the A^^ in ( [4.76|) gives the required A^^ scaling seen in the sum rule. Therefore we 



conclude these additional terms in the shear sum rule are due to the fact operators 
proportional to of dimension 4 in the OPE which acquire non-trivial expectation 
values in the presence of chemical potentials. 



5. Shear sum rule for R-charged M2 and M5-branes 

We will now repeat the analysis and derive the shear sum rules for R-charged M2 
and M5-branes. One of our motivations to do this is to verify if the modifications to 
the sum rules can be explained using the same reasoning employed for the D3-brane 
case. Indeed we will find that that additional terms in the shear sum rule for the 
M2-brane can be explained in terms of expectation values of operators in the OPE 
of the stress tensor. For the M5-brane case we do not find any additional terms in 
the sum rule and this is consistent with the fact that there are no operators of the 
requisite dimension in the supergravity which acquires an expectation value. 
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5.1 M2-branes 

We begin with the analysis of the M2-brane at finite chemical potential. In general 
there are 4 R-charges corresponding to the Cartan's of 5*0(8) therefore it is possible 
to turn on 4 independent chemical potentials. The gravity dual of this system is given 



by the non-extremal R-charged black hole in AdS^ p3[. Using the same methods 
developed for the D3-brane case it is possible to show that the retarded Green's 
function of the T^y component of the stress tensor satisfies property 1. To ensure that 
it satisfies property 2, we need to subtract the divergent terms at large frequencies. 
The constant term which arises in the Green's function as w — t- oo essentially is the 
LHS of the sum rule. Therefore we will directly go ahead and evaluate the divergent 
pieces as w — oo and isolate the term which contributes to the sum rule. 

The metric and the gauge field for the R-charged M2-brane with all the four 
charged turned on is given by 



9u^ \ ti J ju^ 



(5.1) 



4 

Al = -vrTo 
t 3 



2ki JJ(1 + ki) ^ , u = —, Hi = l + kiU 



'H = l[H„ f = n-l[{l + k 

i=l 1=1 

The scalars are given by 

^1/4 

X' = (5.2) 

The four scalars are not independent and constrained by X1X2X3X4 = 1. The above 
background is the solution of the equation of motion of the following action 

S = o.-r. / d x^/^C, (5.3) 

i\2 
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c = R- l{d$r + v{<i>) 



1..^.. .... 1 ' 



2' 4 

where the fields four fields Xi are related to the three independent fields 0j = 

(01,02,03) by 

Xi = exp(-iai ■ 0), (5.4) 
ai = (1,1,1), 02 = (1, -1, -1), 03 = (-1, 1, -l),a4 = (-1, -1, 1). 

The scalar potential is given by 

2 

V{(j)) = —(cosh 01 + cosh 02 + cosh 03). (5.5) 
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Various thermodynamic quantities are summarized as follows: let us first define 

3r_i_ 



To 



47rL2' 

The energy density and pressure are given by 

4 4 



P 



(5.6) 



(5.7) 



i=l 
4 



n 11(1 + h) . 



1=1 



The expression for entropy density and temperature is given by 

3 4 



(5. 



3Wn(l + ^i 



j=i 



The charge density and the chemical potential are given by 

3 



pi 



V2 



7r 



^3/2 rj.2 



\ 



(5.9) 



3 l + k; 



\ 



2ki\{{l + ki). 



i=l 



To obtain the behaviour of the Green's function as a; — ?■ oo, we first write the 
equation of the minimally coupled scalar in the background given by (|5.1| ) in terms 
of Fefferman-Graham coordinates. The metric given in (|5.1| ) admits the following 
form in Fefferman-Graham coordinates. 



1 



-1 + ^ aiz')dt^ + (1 + ^ hiz'){dx^) + dz' 



(5.10) 



i=l 



1=1 



Here we are working with the scaled variables 



iX = — 00. 



(5.11) 



In the Fefferman-Graham expansion given in ( p.lOp it can be shown that the terms 
ai, hi are zero. As we will see subsequently, the asymptotic behaviour of the retarded 
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Green's function is essentially determined by the coefficients 02,^2, 03,^3 which are 
evaluated in (|A.8|) . Now the equation for the minimally coupled scalar in this metric 
( [5 .101 ) is given by 



\ 1=1 / i=l 



0, (5.12) 



where Oj, bi can be obtained given Oj, 6j. For the coefficients of interest, this relation 
is given by 

02 = (a2 - 262), ^2 = «2, 
0-3 = 2(^3 - 263), h = as. (5.13) 

We now rescale the coordinates by defining 

y = \z. (5.14) 

Then the equation for the minimally coupled scalar reduces to 



^ / 00 \ 00 

y \ 1=1 / i=i 



(5.15) 



Note that at A strictly infinity, the equation reduces to that of the minimally coupled 
scalar in the 74^5*4 or the T = background. To solve the equation perturbatively 
in J we define the quantity 

9iy) = ^- (5.16) 
From the expansion of the equation in (|5.15 ), we see that we can expand g as 

11 , , 

9 = 90 + -y^9i + ^5^2 + • • • • (5.17) 

The equations satisfied by these can be obtain from (|5.15|) and the ffist few equations 
are 

2 

^0 + ^0 - -^0 - 1 = 0, 

9'i + 2(fi'o - -)9i - a2y9o - hy"^ = 0, 

y 

g'^ + 2{go - -)g2 - aV^7o - hy^ = 0. (5.18) 

y 

The two independent solutions for the first equation in ( |5.18D are 

^(1) ^ Ki/2iy) ^ y_ ^(2) ^ Ii/2{y) ^ xsmhx 

° Ks/2{y) l + y' ° h/2{y) X cosh X - sinh X ' 
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As in the case of the D3-branes we can solve for the two independent solutions order 
by order in ^ around these solutions. Just as in the case of the D3-brane situation 
we find that the solution around ^'q^'' determines the leading asymptotic behaviour 
of the Green's function as A — )■ cxd. This is because of the condition that at A — )■ oo, 
the solution must reduce to that of the T = case. 13/2(1/) is a growing solution as 
?/ — 7- 00 and thus is not allowed when T = 0. The leading corrections around the 
solution (7q°^ is given by 

2 02 , 2y + 3 



2 '2{l+yy'' 
g, = ^(aa - 2h)y' - ^(^f^(9 + 1% + + (5.20) 
The retarded Green's function is then given by 

hm Gnico, T) = lim ^ + lim - (x'gi'\y) + \g?iy) + gi'\y) 

oj— s>oo oj— >-oo Z4:7l 1j V^O y \ J 

~l~Cr contact ~l~ counter- (5.21) 



Substituting the form for g\ and g\ in the limit we have 

A lim ^ = -\ {a, + 362) , hm ^ = ~{as + 4h). (5.22) 

y^o 2 y^o y^ 8 

From the expressions in ( p.22| ), it is clear that there are two sources of divergence. 
The term proportional to A^ which is identical to the contribution of the T = 
Green's function. The second divergence is due to the term proportional to A. This 
gives a linear divergence in u and we need to subtract this also along with the 
constant term from the limit of g^2^ to ensure that property 2. is satisfied. Therefore 
we define 

contact 

(5.23) 

where 

\/2N^/'^r^ u 1 

= - ^(^2 + 362), (5.24) 

Note that we have subtracted the required terms to ensure that 6Gr{ijJ,T) obeys 
property 2 to apply Cauchy's theorem. The RHS of the sum rule therefore involves 
the following integrand 

lmSGR{uj,T) = p{uj) + 5p{u)- p{uj,0), (5.25) 
= pi^^T) - p{u,0). 
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The LHS of the sum rule is obtained by evaluating 



5GniO,T) = ^^^^^/ -(a3 + &3), (5.26) 

o zio 

Here we have substituted the values of 03, 04 from ( |A.8| ). Thus the sum rule for the 
M2-brane theory in the presence of chemical potential is given by 

"I ^T7 ^(^1 + k2- ks- k4){ki - ^2 + ^3 - ^4)(fci - ^2 - ^3 + ^4) 

o zio 

^ (p»-pT=oM). (5.27) 



Note that p{co) is defined by 

p{u)t = p{oj)t + Sp{u), (5.28) 



where Sp{uj) is given in ( 5.24|) . 



From the above analysis we see that the sum rule is essentially due to the high 
frequency behaviour of the Green's function. High frequency behaviour of the Green's 
function can be extracted from the OPE of the stress tensor. We will now show that 
the linear divergence in u and the additional term involving the chemical potentials 
in the LHS of the sum rule are due to the presence of certain operators in the OPE 
which get expectation values in the presence of the sum rule. The general form for 
the OPE of the stress tensor in this case is given by 

T^.(x)T,.(o) ~ c't^ + i^.p.a/3(x)r„^(o) + 5;,^,(x)a(o). (5.29) 

Here we have included additional terms in the OPE which arise if there are scalars 
denoted by Oa are present in the theory and the three point function(T^yTpo-(9a) is 
non-zero. Note that now the stress tensor is a dimension 3 operator in this case. From 
a similar analysis as in the case of the D3-brane case, we see that a linear divergence 
can be explained due the presence of operators Oof dimension 1 in the OPE. The 
constant terms must arise from operators of dimension 3. Since the expectation 
value of the stress tensor accounts for the term |e in the sum rule, the additional 
terms in the LHS of the sum rule must arise from other operators Oa of dimension 
3. Note that the u divergence in (|5.24|) and the additional term present in the LHS 



of the sum rule in ( |5.27 ) vanish when all the chemical potentials are equal. In such 



a situation from the solution given in (|5.1| ) one sees that there is non-trivial gauge 



field but the scalars given in (|5.2| ) become trivial. Thus the presence of expectation 



values of these operators must be due to the scalars. 
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Let us now examine the scalars in detail. Expanding the scalar potential given 
in ( ^.5|) to quadratic order we find that all the three scalars 0j are of equal mass with 
the mass given by 

m^L^ = -2. (5.30) 

The scaling dimension of the operator dual to the scalars can be found by the mass- 
dimension relation 

A(A - 3) = m^L\ (5.31) 

From this equation we see that the dimension of the operators A can be chosen to 
be 1 or 2. When A is chosen to be 2 the quantization of the scalar is conventionally 
known as the standard quantization. One the other hand when A is chosen to be 
1 the quantization of the scalar obey the alternate quantization. For the neutral 



scalars in the M2-brane background it has been seen earlier ||2^, ^ that one must 
chose the alternate quantization. Let us denote the three operators dual to these bulk 
scalar fields as Oi which have dimension 1 each. We now read out the expectation 
values of these operators in the non-extremal M2-brane background following 
The expansion of 0i close to the boundary is of the form 



a,- + L^P~ + ---. (5.32) 

1^ lyZi 



The factor of L has been re-instated so that /3 has the appropriate dimensions. 
For the alternate quantization, the expectation value of the operators is given by 



the coefficient "^^r-p^ctj- Expanding the scalar fields 0j defined in (|5.4|) near the 



247rL2 

boundary and using ( |5.2|) we obtain 



From these expectation values we see that the candidate operators which are respon- 
sible for the linear uo divergence given in ( ^.24 ) are proportional Of. Indeed taking 



the sum or squares of these operators we obtain 

]{0}) = (5.34) 



2 

i=l i=l 

{{ki + k2-h- kif + (fci - ^2 + ^3 - kif + (fci - ^2 - fcs + kif) 



648 
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In the first line we have assumed large N factorization. Note that the last line in 
( ^■34|) is proportional to the linear omega divergence seen in (|5.24|) . The factors 



of will present in ( |5.34| ) will agree with that present in ( p.24| ) if one assumes 
the structure constant -B^jyp with the normalization of the operators used scales as 
A^~^/^. Indeed this is what is seen from the gravity analysis^. We can also obtain 
the additional term present in the LHS of the sum rule given in ( |5.27D . Note that 
for this we require an operator of dimension 3. We see that the required operator is 
proportional to O1O2O3. Indeed we find 

{0^0202) = {0,){02){02), (5.35) 
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-{ki + k2- k^- ki){ki - k2 + ks- k4){ki - ^2 - ^3 + ^4 



which proportional to the additional term in ( |5.271) with the exception of the N 



scaling. The N scaling can be understood if one assumes that the structure constant 
^liupa scales as A^~^ for these operators. Thus the additional term in the LHS of 
the sum rule ( |5.27| ) is due to these additional operators in the OPE which gain 
expectation value in presence of the chemical potential. 

5.2 M5-branes 

In this section we examine the dual of M5-brane at finite chemical potential and 
finite temperature and re-derive the shear sum rule. Since there are 2 R-charges cor- 
responding to the Cartan's of SO (5) it is possible to turn on 2 independent chemical 
potential,since there are no constrain equation relating this two R-charges. Using 
the differential equation of the massless minimally coupled scalar in this background 
we obtain the retarded Green's function of the T^y component of the stress tensor. 
Examining this differential equation and the same method employed for he D3-brane 
case it can be shown that the regulated Green's function satisfies both property 1 
and property 2 which are necessary for deriving the sum rule. However as will see the 
shear sum rule for the case of M5-brane is not modified. We argue that this must be 
the case since there are no scalars in the M5-brane gravity background which have 
the appropriate dimensions to modify the OPE of the stress tensor. 

The metric and the gauge field for the R-charged M5-brane with all the charges 
turned on is given by p3 



9u \ ti J Ajw^ 



At = -ttTo 
t 2 



\ 



^ 2 
2kiY[{l + ki) ^ , Hi = l + kiu' 
i=i ' 



^ Since we do not have an explicit theory for the M2-branes at present we cannot understand the 
N scahngs from a field theory analysis for the M2-branes. 
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To = |^, e='-^'^N%^llil + t^ (5.36) 



27rL2' 2 

i=l 



2 2 

3 



H, = i + k,u\ n = l[H,, f = n-l[{i + k 



i=l i=l 

The background solution for the two scalars are given by 

^2/5 

X' = (5.37) 

where u = The above background is the solution of the equation of motion of 
the action is 

S = ——- / £xy/^C, (5.38) 



c = R- lid^r + v{<i>) -jYl ^'"'^(p^y 



1..^.. .... 1 ' 



4 

i=l 

Where the two scalar fields X-s are related to = (0i, ^2) b;y 

ai = (v^,y|), al = (-v^,y|). (5.39) 

Note that the two scalars Xi are independent here, unlike the situation in the case 
of the D3-branes and M2-branes. The scalar potential V is given by 

y = -i (^-4X1X2 - 2Xf 1X2-2 - 2Xf 2X3-1 + i(XiX2)-^^ . (5.40) 

The leading terms in the large frequency expansion of the equation for xy component 
of the metric field in the Fefferman- Graham coordinate system is given by 

A" I^p: , 3(a3-5&3)y^ , 

- -(5 + j-6 >^ " (^ + "3^)0 = 0. (5.41) 



where we have defined 



y = \z = A—, iX = —CO. (5.42) 



r r I 



The values of 03, 63 are evaluated in in ( [A. 111 ). Let us define 



03 = 3(03-563), &3 = «3- (5.43) 
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Since the Green's function is proportional to g = ^, we can write the equation in 
( CT) as 

9' + 9'-k5 + - (1 + &3^) = 0. (5.44) 

We then find the solution to the leading order by writing g as 

9 = 90 + ^9i + ■■■■ (5.45) 

The equations satisfied by each of the terms in the expansion of g can be obtain from 
( p. 44 ). They are given by 

5 
— ( 



fi'o + 5-0 - -fi'o - 1 = 0, 



5 _ 

g[ + 2gogi - -gi - ~a^y^go - hy"^ = 0. (5.46) 

y 



The relevant solution for our analysis are 



^° " KsivY " 2 + 14 Kt ^^-^^^ 

The retarded Green's function is then given by 

hm Gniu,T) = lim lim (x^g^'\y) + g['\y) 

+ G'contact + Ccounter- (5.48) 

The term proportional to which is identical to the contribution of the T = 
retarded Green's function. It is divergent in the limit A — )■ oo. To ensure that 
the property 2 is satisfied we need to subtract the divergent piece and the constant 
pieces from the Green's function to write the regulated Green's function. Therefore 
we define 

6Gr{u, T) = Gr{u, T) - Gr{u, 0) - G 

contact 

(5.49) 

^3 15^6 

+ 3^14Z^^^^^^ + ^^)- 
The LHS of the sum rule is obtain by evaluating 

= ^e. (5.50) 

Here we have substituted the values of 03, 63 from (|A.11| ). Thus the sum rule for the 
M5-brane in the presence of chemical potential is given by 



3 

r 



-[ ^ (p(uj) - pT=o{i^)) . (5.51) 



-40- 



Note that for the case of R-charged M5-branes there is no additional contribution 
in the shear sum rule compared to uncharged case. This fact can be understood as 
follows. Let us examine the masses of the scalars in the theory. From expanding the 
potential for the scalars in ( ^.40|) we see that the two scalars (pi in the theory have 



mass given by by 

m^L^ = -8. (5.52) 

From the mass-dimension relation we see that there are two possible choices of mass 
for the operators dual to these. They are A = 2 or A = 4. The stress tensor for the 
6-dimensional theory of the M5-brane is an operator of dimension 6. Thus if there 
are terms which are finite in the large frequency expansion of the OPE, they must 
be of dimension 6. The general form for the asymptotic expansion of the scalar field 
in this case is given by 

<^^^ = ^ + -, + ■■■■ (5.53) 
From examining the background value for the scalar field given in (|5.37| ) we see that 



for both the scalar fields, = and /Sj 7^ 0. In the standard quantization of the 
scalar fields, the dual operators have dimension A = 4 and the expectation value of 
the scalar field is proportional to the value of But in this situation there is no 
operator which has a dimension A = 6 which gains expectation value. This implies 
there is no finite term in the large frequency expansion of the OPE of the stress 
tensor ^. If we choose the alternate quantization, then the operators dual to the 
fields (pi have dimension 2, but from the expansion of the scalar field we see that 
their expectation value is proportional to the value of which is zero. Thus in the 
alternate quantization the operator dual to the scalar field does not gain expectation 
value. Therefore it cannot contribute to either a divergence or the finite term in the 
OPE of the stress tensor. These arguments imply that there is there is no correction 
to the sum rule which is consistent with our explicit calculation. 



6. Conclusions 



By examining the shear sum rule for the case of R-charged D3-branes and M2- 
branes we arrive at the general observation that sum rule acquires additional terms 
proportional to the chemical potential. These additional terms are the result of 
expectation values of operators of dimension 4 and dimension 3 for the case of D3- 
branes and M2-branes respectively. For the case of the M5-brane this no change in 
the sum rule at finite chemical potential. This is explained by the fact that there is 
no operator of dimension 6 whose expectation value is turned on at finite chemical 

^The operator dual to the scalar field with dimension A = 4 can give rise to a divergence if 
there the 3 point function of this operator with the stress tensor {T^^Tp^O) is not zero. From our 
gravity analysis since there is no divergence it is clear that this three point function vanishes. 
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potential. Our analysis indicates that the LHS of the sum rule contains information 
of the OPE coefficients of the stress tensor with scalar operators of the appropriate 
dimensions. Thus this approach can be used to compute three point of functions of 
the stress tensor with these operators from gravity. This will involve writing down the 
relating the coefficient B'^^p^ given in ( |4.60| ) in terms of the coefficient which occurs 
in the three point function using a similar analysis of []T^. It will be interesting to 
explicitly check the three point function obtained in this manner against a direct 
evaluation three point functions from gravity. 

Another direction to extend this work would be to examine sum rules for the 
R-charge correlators in = 4 Yang-Mills at finite chemical potential. It will be 
interesting to see if the effects of triangle anomalies can be seen from the sum rule 
just as hydrodynamics was modified on the account of triangle anomalies ^6). Finally 
we mention that chemical potential is an important parameter in the QCD phase 
diagram and it will be interesting to see if sum rules can be derived in this context. 
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A. Fefferman- Graham coordinates 

In this appendix we will detail the co-ordinate transformation which takes the metrics 
discussed in this paper to the Fefferman- Graham coordinates. We will see that for 
the non-extremal R-charged D3 and M5-branes, the metric in Fefferman- Graham 
co-ordinates takes the form 



-. / oo oo 

ds' = -;^l (-1 + J2<^,{z"y)dt' + (1 + J2H^'"y)id^') + dz" 1 . (A.l) 

\ j=l i=l / 



For the case of the non-extremal R-charged M2-brane the metric in Fefferman- 
Graham co-ordinates takes the form 



-. / oo oo \ 

^ (-1 + 5^a,(z')0rft^ + (1 + Y.^^{z'r){dx') + dz'^ , (A.2) 
^ ^ \ i=i i=i J 



here L = L/r+. In this appendix we will evaluate the coefficients ai,bi which are 
necessary to obtain the sum rule. 
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Non-extremal R-charged D3-brane 

The metric R-charged D3-brane metric given in in ([4.1|) can be transformed to the 
Fefferman- Graham form by first redefining the co-ordinate z = ^ in terms of the 
Feffer man- Graham co-ordinate z' as a power series 

z = Q{z'), (A.3) 
= z'{l + az'^ + /3z'* + 0{z'^)). 

After expanding the coefficients of the metric ( [4.1| ) as a power series in z around 
z = we substitute the expansion ( |A.3| ) in the metric ( |4.1j ). We then determine the 
coefficients a, /3, ■ ■ ■ by requiring that the coefficient of the dz'^ term reduces to 
Performing this procedure to O(z^) results in the following values of a, P 

a = l{h + k2 + h), (A.4) 
o 

/3 = + kl + kj) - {kik2 + k2k^ + k^ki) - 9(1 + ki + k2 + k^ + kik2k^)}. 

Finally we substitute these values back in the remaining coefficients of the metric to 
obtain 

Qi = 0, hi = 0, (A.5) 
as = 7^{27(1 + ki + k2 + k^ + kik2k3) + 25{kik2 + ^2^:3 + k^ki) 

+ 2{kl + kl + kl)}, 
b2 = 7^{9(1 + ki + k2 + k3 + kik2ks) + ll{kik2 + ^2/^3 + k^ki) 

-2{kl + kl + kl)}. 
Non-extremal R-charged M2-brane 

We now perform the same procedure for the metric given in (|5.1|) . Let us define the 
co-ordinate z as 

z = n{z'), (A.6) 
= z'{l + az' + (3z'^ + 7z'3 + 0{z'^)). 

Demanding that the coefficient of dz'"^ to be of the form of ^ to 0{z'^) in the new 
coordinate system determines a,/?, 7. They are given by 

a = ^{ki + k2 + k^ + ki), (A.7) 
/3 = hiikl + kl + kl + kl) + 10(/ci/c2 + kik^ + k^k^ + k2k^ + ^2^:4 + k^ki)], 

i = ^^m+ki+ki+ki) 

+ ll{klk2 + klki + klkj, + klki + k^k^ + klki + kjk^ + klk2 + klk^ + klk2 + klk^ + kjks) 
- 2{kik2k3 + kik2ki + kik^k^ + /c2^3fc4) - 64(1 + ki + k2 + k^ + k^ 
+ kik2 + kiks + klki + ^^2^3 + ^2^4 + ^3^4 + fciA;2fc3^4)}- 
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We then substitute these values back in the other coefficients of the metric to deter- 
mine the values of a^, bi which results in the following values. 

a^ = bi= 0, (A.8) 
«2 = ^{^{kl + kl + kl + kl) - 2{kik2 + k^k^ + kik^ + fcz^s + k2k4 + k^k^)}, 

as = ^{-(^1 + kl + kl + kl) + lA{kik2h + kik2ki + kik^k^ + k2kj,ki) 

+ {klk2 + klki + klk^ + klki + klk^ + klki + klk^ + klk2 + klk^ + klk2 + A;^A;4 + fc^Ajg) 
+ 16(1 + ki + k2 + k'i + k4 + kik2 + kik^ + kik^ + ^2/1:3 + ^2^4 + ^3^4 + ^1/^2 ^3/^:4)}, 

62 = -^{3(A;i + fc^ + A;^ + kl) - 2{kik2 + fci/cg + kik^ + ^2^:3 + ^2^4 + k^ki)], 

63 = ^{(^? + A;^ + + A;4) + 10(A;iA;2A;3 + /i;i/;:2/i;4 + /i;i/^3/i;4 + k2kj,k^ 

- (A;^A;2 + k\ki + A;^A;3 + A;^fci + k\ki + + klkj, + A;^A;2 + A;^A;4 + klk2 + fc^A;4 + klk-^) 
+ 8(1 + fci + A;2 + /i;3 + A;4 + A:i/i;2 + k^k^^ + A;iA;4 + ^2/^3 + A;2A;4 + ^3^:4 + k\k2k'iki^'] . 

Non-extremal R-charged M5-brane 

We follow the same procedure again. We ffist expand the metric around z = 0, 
retaining terms to order Oi^z^) and then apply the coordinate transformation 

z = ^{z'), 

= z'{l + az'^ + (3z'^ + -iz^ + 0{z")). (A.9) 

By demanding that the coefficient of dz'"^ is of the form of ^ to 0{z''^) in the new 
coordinate system we can determine and 7 . They are given by 

„ . ki + k2 (1 + A;i)(l + A:2) 
a = 0, (5= , 7 = • (A.IO) 

Substituting these values back in the metric we determine the other coefficients of 
the metric to be given by 

ai = hi = 0, 02 = 62 = 0, (A. 11) 

as = ^(1 + ki){l + k2), h = 1(1 + h)il + k2). 
o 
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